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CHAPTER 1 


INTRODUCTION 

For several years the NASA Langley Research Center has been flying 
an instrumented F106B aircraft into thunderstorms collecting information 
on direct lightning strikes. The ultimate goal of the effort is to be able 
to characterize the lightning environment so that the effect of direct 
lightning strikes to arbitrary aircraft can be predicted. The F106B carries 
instruments to measure such quantities as electric and magnetic fields on 
the aircraft, current flowing on the pitot boom, and ambient temperature. 

The electromagnetic measurements are made as a function of time with a 
resolution of 10 nanoseconds. This is sufficient to resolve the evolution 
of a lightning channel as it forms, usually near the extremities of the air- 
craft. The task of the analyst is to use the electromagnetic responses measured 
on the F106B to determine the characteristics of the lightning which caused 
the response. 

This report concentrates on the 1982 direct strike data [1]. The ana- 
lytic thrust is the investigation of triggered lightning on the F106B, because 
recent studies [2] have shown this to be quite likely. The major analytical 
tool is the computer code T3DFD, a three dimensional finite difference model 
which solves Maxwell's equations in the time domain. The physics of the air 
is taken into account with a three species air chemistry model which provides 
a nonlinear air conductivity to Maxwell's equations. 

A block model of the F106B, which is how the airplane appears to 
the finite difference code, is shown in Figure 1.1. The locations of the 
B-dot and D-dot electromagnetic sensors are indicated in the figure. The grid 
resolution of the finite difference model is one meter in the direction along 
the fuselage and one-half meter in the other two dimensions. This requires a 
time step of one nanosecond, which is the temporal resolution of the model. 

The specifics of the finite difference model have been covered in some detail 
in previous reports [3,4], and will therefore be omitted here. 

As stated earlier, the body of this report is concerned with the 
interpretation and analysis of the 1982 direct strike lightning data in light 
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Figure 1.1 Block Model of F106B Showing Sensor Locations 


of the possibility that much or all of it may be from triggered lightning. 
Chapter 2 is concerned with data classification in the time domain, searching 
for trends that could indicate similar physical processes. Chapter 3 analyzes 
the data in the frequency domain, with emphasis on aircraft resonances which 
should be present. Chapter 4 deals in a general fashion with the triggered 
lightning environment, and its relation to the direct strike program. Chapter 5 
discusses the modeling of triggered lightning on the F106B and presents results 
from the computer code T3DFD. Chapter 6 presents a new enhanced nonlinear 
model which includes fluid equations for particle momentum and energy densities. 
This model, because of its size, is currently only two dimensional and is still 
under development. Chapter 7 is concerned with the details surrounding sub- 
grid implementation. This allows one to model a given region in a finite dif- 
ference code with greater spatial resolution than is present in the rest of 
the problem space. Chapter 8 details the calculation of transfer functions for 
the F106B in three different current injection configurations. Chapter 9 pre- 
sents a summary and conclusions. 
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CHAPTER 2 


DATA CLASSIFICATION 

In an effort to identify trends in the measured data from 1982, 
the D-dot and B-dot records have been classified into categories according 
to waveshape. The classifications are by necessity somewhat subjective, and 
it is certainly possible that some waveforms belong in a different category. 
The categories are intended to highlight waveforms with similar character- 
istics, which may indicate similar physical phenomena. Descriptions of the 
eight D-dot categories and six B-dot categories are given below, with the 
number of waveforms in each category in parentheses. 

D-Dot 

1. Slow positive, then negative, bipolar pulse. (35) 

2. Double bipolar pulses separated by ^200-400 ns, initially >0, with 
f D dt i$. (11 ) 

3. Single negative pulse with little structure. (15) 

4. Bipolar pulses, first negative, then positive (3) 

5. Fast bipolar pulses, initially positive. (21) 

6. Initial large fast pulse with much structure and having a late time tail, 
mainly unipolar. (24) 

7. Simple unstructured positive pulse. (6) 

8. Miscellaneous. (8) 

B-Dot 

1. Initially positive pulse with much late time structure. (7) 

2. Slow positive and negative half cycles, usually small. (1) 

3. Initially positive single cycle followed by low amplitude structure. (15) 

4. Unstructured single positive cycle. (5) 

5. Single negative slow pulse. (3) 

6. Fast double pulses, separated by 200-500 ns, initially positive. (19) 

The waveforms themselves along with their Fourier transforms are 
presented in Appendix A. The format of the majority of the figures there is 
to overlay three of the time domain waveforms in a given category in the 
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upper left corner of the page. The legend giving the flight number and 
characteristics of the strike (attach point, temperature, flight altitude) 
is coded according to the line type of the waveform (dashed, dotted, or 
solid). The other three positions on the page are the Fourier transform 
of the time domain waveforms coded as in the time domain plot. The exceptions 
to this format occur if there are not three time domain waveforms to over- 
lay, as in Category 8 of the D-dot records (Miscellaneous). The format of 
the first exception is transparent. For Category 8 of the D-dots, no over- 
lays were made because the waveforms were judged to be unique. Hence the 
two plots on the left hand side of the page are time domain waveforms with 
the corresponding Fourier transforms on the right of the page. 

Not all of the 1982 data records are included in this cassifi- 
cation. Records having very small amplitude (vtwo or three digitization 
increments) and records with saturated sensors were omitted. Although the 
saturated records probably could have been included in the time domain 
classification, it was impossible to get accurate Fourier transforms of them. 

Finding a correspondence between the classification categories 
and particular physical phenomena is difficult. There appears to be no ob- 
vious correlation between temperature and category or flight altitude and 
category. The correlation between attach point and category may be slightly 
better, although the attach point for the majority of strikes is unknown. 

It is very possible that D-dot records from Categories 2 and 5 and B-dot 
records from Categories 3 and 6 are caused by lightning triggered at the 
nose of the aircraft. (See Chapter 5 for representative waveforms). 

These categories probably represent the same type of phenomena, except 
that in the categories with double pulses, a similar event happens twice 
(e.g. air breakdown at nose and later at tail). The hypothesis that triggered 
lightning causes these pulses is supported by the simultaneous D-dot and 
B-dot records. The models for natural and nearby lightning do not reproduce 
the relative amplitude and timing of D-dot and B-dot responses seen in the 
measured records, while the model for triggered lightning does. Although it 
is not possible to exclude natural and nearby lightning as responsible for 
the measured records without examining in detail all possible cases, the 
evidence so far tends to lend credence to the triggered lightning hypothesis. 
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The other categories have not been identified with physical events 
as yet, although if triggered lightning is responsible, it could be triggered 
at points other than the nose of the aircraft. Triggering at locations other 
than the nose will be investigated in the near future and may give insight 
into the physical events which produce many of the D-dot and B-dot categories. 
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CHAPTER 3 

STATISTICAL ANALYSIS OF FOURIER TRANSFORMED DATA 
3.1 Fourier Transform Analysis of Measured Data 

From the time domain data measured on the F106B it is difficult 
to recognize the presence of aircraft resonances. This can be much more 
easily done if the data is transformed into the frequency domain. This has 

been accomplished for most of the measured data, and a statistical study 
performed on the results. 

The data transformed was that from the D-dot forward, B-dot longi- 
tudinal, I (current on nose boom), and I-dot (time derivative of current 
on nose boom) sensors. The only records from these sensors which were not 
transformed were records with very small amplitude (* few digitization 
steps), obviously saturated records, and records having data dropouts. 

The magnitude of the transform of all the B-dot and D-dot records are 
presented in Appendix A, along with their associated time domain waveforms. 

I and I-dot records are not included in Appendix A because their 
numbers were not sufficient to warrant a classification as was done with the 
D-dot and B-dot records. 

The statistical analysis of the Fourier transforms proceeded as 
follows. For each transform magnitude a search was done to find the maxima 
and minima of the transform as a function of frequency. A database was 
constructed which included the frequency location of all the maxima and 
minima as a function of type of record (D-dot, B-dot, I, I-dot). The database 
also included the absolute value of the transform at the maxima and minima. 
Because the absolute value of a transform at any point is not a particularly 
useful parameter in searching for resonances, another parameter was de- 
veloped which measured the relative importance of a given maximum or mini- 
mum in a frequency spectrum. If one calls the maxima and minima of a spectrum 
the resonances , this parameter will be called the "relative amplitude" of 
the resonances. For the definition of the relative amplitude consider 
Figure 3.1. In that figure is shown a typical piece of a Fourier transform, 
with frequency on the horizontal axis and the magnitude of the transform 
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on the vertical axis. The relative amplitude of the peak at frequency f^ 
is defined to be the average of the height of the peak above the two minima 
on either side. 

Relative amplitude of peak at f. is given by: 

1 (3.1) 

Relative amplitude = 1/2 (Tj-Tg+T-j-Tj). 

Similarily the relative amplitude of a minimum is defined to be its 
average depth below the two peaks on either side of it. Using these 
definitions the relative amplitude parameter is always a positive number 
and gives a measure of the prominence of a resonance in a given spectrum. 

Having defined the relative amplitude parameter, it was then in- 
corporated into the database and used as a filter in searching for air- 
craft resonances. This was done by producing histograms showing number of 
resonances in a given frequency bandwidth as a function of frequency, 
type of record, type of resonance (maximum or minimum), and relative ampli- 
tude. As mentioned the relative amplitude was used as a filter. First all 
resonances of a given type were included in a histogram. Next only those 
with relative amplitudes greater than 2.5 dB were included. The relative 
amplitude filter was raised in this way until only a statistically insig- 
nificant number of resonances could meet the criteria. Using this method 
it was hoped that as the relative amplitude filter was increased, physical 
aircraft resonances would rise out of the noise in the spectrum. The re- 
sulting histograms are presented in Figures 3.2 to 3.39. I-dot sensor results 
are not presented, because the small number of I-dot records was statistically 
insignificant. The histograms have a bandwidth of 500 kHz, so that each 
bar represents the number of resonances falling In a 500 kHz window centered 
on the frequency label below it. For example, the bar labelled 15 MHz 
actually contains all the resonances falling between 14.75 MHz and 15.25 MHz. 

Interpretation of the histogram figures is relatively easy. As 
one moves from small relative amplitude filter to larger relative amplitude 
filter, the total number of resonances goes down. This occurs as the small 
resonances produced by numerical and measurement noise are eliminated. There- 
fore the signal to noise ratio of the remaining data is improved. Using the 
D-dot maximum histograms as an example, one starts with a relatively flat 
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figure 3.35 i-dot maxima, amplitude filter= 
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graph for a zero dB filter. By the time one reaches the 10 dB filter graph, 
several frequencies have begun to assume prominence. There is a low fre- 
quency bulge below 5 MHz, another between 7 and 10 MHz, and indications of 
resonant behavior in the ranges around 17 MHz, 21 MHz, and 27 MHz. At even 
larger relative amplitude filters the 7 to 10 MHz and 21 MHz resonances 
become dominant. It may be noted that the period of the 7 MHz oscillation 
(% 140 ns) corresponds to the round trip travel time of a signal from the 
nose of the aircraft to the tail and back. 

The B-dot histograms are less amenable to simple explanations, 
possibly because of the location of the sensor. At the largest relative 
amplitude filter values, there does appear to be some activity near 7 MHz, 
but it does not rise above the background as prominently as in the D-dot 
histograms. This may be because the D-dot sensor, located forward on the 
fuselage, is affected most strongly by fuselage resonances. Although the 
B-dot sensor can also be affected by these, it will also respond to modes 
which include the wings, and then may be expected to contain frequencies 
other than the basic fuselage resonance. 

If this is the case, one would expect the current histograms to 
exhibit similar behavior to the D-dot' s. This is indeed the case, as one 
sees the 7 to 10 MHz resonance rising above the background for the largest 
filter values. Also the resonance near 21 MHz is prominent again in the 

current data. 

This method of analyzing the Fourier transformed in-flight data 
has shown that aircraft resonances do indeed exist in the data, and that 
they are predominantly fuselage resonances. This is consistent with a 
source at the nose of the aircraft, which is, of course, where many of the 
strikes are known to have attached. 
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CHAPTER 4 


TRIGGERED LIGHTNING CONCEPTS 

4.1 Introduction 

Triggered lightning is so-called because it occurs only in the 
presence of a conducting object which is able to enhance the local electric 
field enough to cause air breakdown. Although any conducting object will 
produce some field enhancement, the largest enhancements will occur for 
objects having sharp points or edges, particularly if those points and edges 
are oriented in the direction of the ambient field. An aircraft in flight 
such as the F106B has many sharp points (e.g. nose) and sharp edges (e.g. 
wings, tail) around which the fields will be enhanced if the aircraft is 
immersed in a static electric field. These locally (within a half meter or so 
from the aircraft) enhanced fields are likely to be a factor of ten or more 
larger than the ambient field. Hence it is considerably more likely that 
initial air breakdown and formation of a lightning channel will occur in the 
presence of the F106B than in its absence. 

Since the possibility exists of triggered lightning occurring during 
the F106B flights, several questions must be addressed. First, what is the 
relative likelihood of a given strike being triggered or natural lightning? 

Here natural refers to the event in which the F106B accidentally happens to 
be in the path of a lightning stroke which would exist without the aircraft. 
Second, if one assumes that triggered lightning is important, what can be said 
about its variation with flight parameters such as altitude and flight heading? 
A related issue concerns how one could modify the flight profile to increase 
strikes, especially at low altitudes. Third, the inflight data must be eval- 
uated with the possibility that much of the data is from triggered strikes. 
Finally, the triggered lightning hypothesis must be applied to the over- 
all goals of the program. That is, how can triggered lightning data gleaned 
from the F106B be used to determine the effects of lightning, triggered 
or natural, on other aircraft. 

4.2 Relative Likelihood of Triggered and Natural Lightning 

To determine the relative importance of triggered versus natural 
lightning, it is helpful to investigate the circumstances under which both 
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occur. The requirements are simply stated for triggered lightning. The static 
electric field which the aircraft sees must be large enough and oriented 
properly so that the locally enhanced fields at the extremes of the aircraft 
exceed the local air breakdown value. The ambient field must also be large 
enough to form and propagate a channel , once a breakdown region has been 
established. A net charge on the aircraft can either suppress or enhance 
the breakdown, and this may significantly change the requirements for the 
occurrence of triggered lightning. Also a net charge can change the location 
of initial air breakdown by altering the local field distribution around the 
aircraft. Finally, one may conclude that triggered lightning will occur if 
the proper conditions are satisfied. Essentially what this means is that the 
local field somewhere on the aircraft must exceed the breakdown values of 
approximately 1.5 x 10 6 V/m at 20,000 ft and 3 x 10 6 V/m at sea level. It 
should also be noted at this point that triggered lightning will occur at 
the lowest possible field level. Hence all the strikes should be similar 
in amplitude depending on where on the aircraft the initial breakdown occurs. 

The conditions under which natural lightning will interact with an 
aircraft are much less certain. The element of chance is far more important, 
because the aircraft must at least be "near" a normal lightning channel in 
order for attachment to occur. The term "near" requires some definition. 

The aircraft must be close enough to perturb the electric field distribution 
of the channel, to alter the normal channel path and cause it to pass through 
the aircraft. A general rule of thumb for the perturbing effects of a con- 
ducting object is that they extend the maximum dimension of the object in 
all directions. That is, an aircraft can perturb the fields in a volume about 
three aircraft lengths in diameter. The region of strong perturbation is much 
smaller than that. For example, the local field around the pitot boom of the 
F106B can be quite large, but its extent is very small, reaching only a meter 
or two from the boom. Beyond that the field is essentially the ambient field 
slightly enhanced by the presence of the aircraft as a whole. Therefore the 
large locally enhanced fields around the extremes of an aircraft are likely 
to have little effect on natural lightning channels. This means that natural 
lightning strikes to the nose (and other sharp points and edges) are only 
slightly more probable than strikes to other parts of the aircraft. Then in 
conclusion one can say that natural lightning strikes to the F106B will occur 
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only if the channel happens to intersect a small volume around the aircraft. 
In addition the statistics of attachment may not be strongly biased toward 
edges, (i.e. nose, tail, wing tips), but may be more evenly distributed 
around the aircraft. 


A third type of strike to the aircraft may be possible, which will 
be called a "hybrid" strike. In this case, a natural leader channel would be 
positioned somewhere near the aircraft, but not near enough for the aircraft's 
field distortion to alter the channel's course. The lightning leader leaves 
behind it in the channel a certain charge per unit length, which produces 
to lowest order a static field around the channel. This static field, at the 
location of the aircraft, may be large enough to cause a triggered streamer 
from the aircraft. The streamer would then certainly propagate to the natural 
channel, becoming a branch of the natural stroke. It is unclear whether a 
strike of this kind would be more characteristic of a natural strike or a 
triggered strike. 


A simple calculation can be done to get some idea of how near a 
natural stroke must be to the aircraft to cause a hybrid strike. For simpli- 
city one may assume that the natural channel is vertical and has a uniform 
charge density per unit length. Because only order of magnitude numbers are 
desired, it will also be assumed that the situation is quasi-steady state. 
That is, changes in the field distribution around the channel occur slowly 
and the static solution is approximately correct. The electric field around 
a line charge of charge per unit length X, which is effectively infinitely 
long, is radially directed and given by: 


E = — - — 
r 2tte r 


(4.1) 


Choosing an enhancement factor for the aircraft of about 10, which 
is typical around the nose of the aircraft, one finds that triggering will 
occur at a distance. 


r trig 




(4.2) 


where represents the air breakdown field. A typical value for x, calcu- 
lated for a channel lowering 5 coulombs of total charge over a 3 kilometer 

-3 

distance, is 1.667 x 10 coul/m. E. varies with altitude, at sea level 

c b ’ 

being about 3 x 10° V/m and at 20,000 feet about 1.5 x 10° V/m. Substituting 
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these numbers into (4.2) gives triggering distances of approximately 100 
meters at sea level and twice that at 20,000 feet altitude. This is a sig- 
nificantly greater distance than the distance at which the presence of the 
aircraft can directly alter the course of a natural lightning channel. Hence 
this hybrid strike is much more probable than a natural strike. It may even 
be that the natural strike as defined here does not exist in that all strikes 
involving a natural lightning channel occur according to the hybrid scenario. 
If this is the case then natural strike attachment locations will be strongly 
biased toward sharp points and edges, as these are the locations from which 
hybrid streamers will emerge. 

It is possible that the hybrid process could account for some of 
the variability in the measured data. The reasoning behind this is the fol- 
lowing. The normal triggered lightning model which has the aircraft flying 
into a slowly increasing ambient electric field requires that the triggering 
occur at the smallest possible field level. This must result in aircraft 
responses which are very similar in amplitude and structure, because the 
only really significant adjustable parameters are the orientation of the 
field with respect to the aircraft and the net charge on the aircraft. How- 
ever, in the hybrid model the natural channel can appear at any distance 
from the aircraft in a relatively short time, so triggering may occur at 
field levels which are significantly above the minimum level. This may re- 
sult in aircraft responses which have a much wider variety, both in amplitude 
and structure. 

Given the conditions under which natural and triggered lightning 
strikes occur, it appears far more likely that most strikes to the F106B 
are triggered. Triggered lightning is certain to occur under the proper 
conditions, while natural lightning is probabal istic even under ideal con- 
ditions. 


4.3 Triggered Lightning Environment 

To investigate the triggered lightning environment a typical 
thunderstorm will be used [5], The static electrical characteristics of 
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this thunderstorm are a +40 Coulomb charge centered at 10 km above ground, 
a -40 Coulomb charge centered at 5 km, and a +10 Coulomb charge centered 
at 2 km, as shown in Figure 4.1. The model of Figure 4.1 will be simplified 
slightly by assuming the charges can be considered to be point charges, 
and that they are located in a vertical line. For this simplification the 
static electric field can be calculated as a function of space. The cal- 
culation assumes that the earth is a perfect ground plane, so it can be 



Figure 4.1 Typical Thunderstorm Charge Distribution [5]. 

P = + 40 coul., N = - 40 coul . , p = + 10 coul. 


replaced by image charges. The results of the calculation are shown in 
Figures 4.2a - 4.2c. These figures show contours of constant field as a 
function of space. The vertical scale is altitude and the horizontal scale 
is radial distance from the (assumed) vertical line of charges. Figure 4.2a 
shows contours of the radial component of the field. Figure 4.2b the verti- 
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Contours Showing Variation of the Radial Horizontal 
Electric Field Around a Typical Thunderstorm Charge 
Distribution. Contour Labels Indicate Electric Field 
Strength. The Dashed Lines Indicate th f Ration of 
the Minimum Field Necessary to Trigger Lightning on 
the F106B. 
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Figure 4.2b Contours Showing Variation of the Vertical Electric 

Field Around a Typical Thunderstorm Charge Distribution 
Contour Labels Indicate Electric Field Strength. The 
Dashed Lines Indicate the Location of the Minimum Field 
Necessary to Trigger Lightning on the F106B 




Figure 4.2c 


Contours Showing Variation of the Total Electric 
Field Around a Typical Thunderstorm Charge 
Distribution. Contour Labels Indicate Electric 
Field Strength 
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cal component of the field, and Figure 4.2c the total field magnitude. It 
should be remembered that the contours shown on the figures were calculated 
assuming point charges, so the fields shown are really upper limits on the 
actual field. If one of the contours is inside an extended charge cloud, 
the actual field is likely to be somewhat smaller at that point than the 
contour indicates. 

Figure 4.2a is the most interesting for the purposes of investi- 
gating triggered lightning on the F106B. Assuming the aircraft is always 
in essentially level flight, the horizontal field is the one which will be 
most enhanced. Linear calculations done on the F106B show that fields 
along the axis of the aircraft fuselage are enhanced by about a factor 
of ten at the nose. Wing - wing ambient fields are enhanced by about a 
factor of seven at each wing tip. Vertical ambient fields are enhanced by 
about a factor of three at the tip of the vertical stabilizer. Hence for 
level flight vertical ambient fields are much less important than radial 
fields, unless the ambient vertical field is several times larger than 
the radial field. 

Another factor that must be considered is that the minimum air 
breakdown field is larger at lower altitudes because of the increased air 
density there. Hence a field which causes triggered lightning at a high 
altitude (% 10 km) may not do so at lower altitudes because of the greater 
air density. The dashed lines in Figures 4.2a and 4.2b illustrate this 
effect. They are the locus of points in space at which an electric field 
of breakdown strength can be reached on the aircraft with proper orientation. 
They also indicate the regions inside which triggered lightning can occur 
for the F106B in level flight. Note that the volumes are much larger at 
high altitude, extending to almost two kilometers from the charge center 
for the radial field case. At the very low altitudes ('v, 2 km) the trigger 
region is only a few hundred meters across. In fact because the charges are 
really extended and not point sources, it is possible that at these low 
altitudes no triggering region exists at all. 

This very simple model illustrates the environment for triggered 
lightning in the vicinity of a typical thunderstorm. Many complicating factors 
have been left out. In a real thunderstorm the charge centers are not points 
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but exist in extended regions. They also may be clumped into several smaller 
charge centers, around which locally high electric fields may exist. Second- 
ly, in most thunderstorms there is a shear in the vertical direction, so 
that the charge centers are not in a vertical line. This destroys the 
cylindrical symmetry. Third, other thunderclouds in the vicinity will 
alter the field distribution around any given cloud. Finally, the effect 
of a net charge on an aircraft has been omitted. Therefore, the preceding 
discussion should be regarded as a guide only to the generally expected 
behavior of triggered lightning and not to any specific thunderstorm 
situation. 

4.4 Triggered Lightning and the Measured Data 

In interpreting F106B direct strike data, it is of some interest 
to investigate in a simple way the expected response of the aircraft. That 
is, given that triggered lightning occurs, what will the measured D-dot 
and B-dot records show, and can one expect any significant difference from 
records produced by natural lightning ? 

To begin, consider the situation in Figure 4.3. This would correspond 
to the F106B flying directly toward a positive charge center. The aircraft is 
polarized oppositely to the ambient field, resulting in a field at the D-dot 
sensor point which points toward the fuselage. Presumably this field has 
grown slowly over several seconds, as the aircraft flies into a slowly in- 
creasing ambient field. This slow growth is far below tne trigger level for 
the derivative sensors on the aircraft. For this geometry the largest enhanced 
field is at the nose of the aircraft, and eventually the enhanced field there 
will become large enough to cause corona. At that time electrons will flow 
off the nose causing the field at the D-dot sensor location to become less 
negative. This will happen very quickly, and will probably trigger the on 
board instrumentation. As the air breaks down at the front of the F106B, the 
field will increase at the aft end, possibly causing another breakdown there. 
That event is dependent on the characteristics of the breakdown at the nose. 

The corona at the nose effectively increases the length of the F106B and 
therefore enhances the local field at the rear of the aircraft. Therefore a 
large corona at the nose will be more likely to cause a breakdown behind 
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Figure 4.3 Simplified Drawing of the Expected Behavior of 
Triggered Lightning on the F106B for Ambient 
Electric Field Oriented Nose to Tail 



the plane also. It should also be mentioned that the hot plume of exhaust 
behind the aircraft may also increase the probability of an air breakdown 
there. This has yet to be investigated in detail. 

Figure 4.3 also shows the expected behavior of the fields and 
their time derivatives at the D-dot and B-dot sensor positions. The plots 

are simplified to the extent of ignoring any behavior caused by the re- 
sonances of the aircraft. Actual records could be expected to look like 
those in Figure 4.3 only with respect to general features. The pulses seeen 
there could be separated by several hundred nanoseconds or longer, depending 
on the growth of the corona at the nose. It is significant that both single 
pulses and double pulses as shown in Figure 4.3 are seen in the measured 
data. (Consider D-dot categories 2 and 5 and B-dot categories 3 and 6 of 
Appendix A). 

An example from the measured data that may illustrate this behavior 
is shown in Figure 4.4. Presented there are the D-dot and B-dot records from 
Flight 82-041 run 20. There are two D-dot and B-dot events for that record 
separated by about 170 microseconds, which behave in a similar manner to 
the simplified waveforms shown in Figure 4.3. This is good evidence that 
the events depicted in Figure 4.3 are actually occurring during the flights 

of the F106B. 

Figure 4.5 shows the D-dot and B-dot records from another event 
(82-041, run 6) in which some different process is occurring. Here, the 
records happen about 700 microseconds apart, and their general behavior is 
different than that of Figure 4.3. This type of measured data, in which two 
or more D-dot and B-dot events exist on a single data record, could be 
quite helpful in characterizing the lightning phenomena occurring on a given 
flight. For example, two D-dot events on a single record could indicate 
breakdown at two different locations on the aircraft. The relative amplitude 
of the two events can provide information as to the location of these break- 
downs. From that information, in turn, can be inferred possible ambient field 
strength, orientation, and aircraft charge. 

The case in which the ambient field is reversed in direction is 
shown in Figure 4.6. This would occur if the aircraft were flying away 
from a positive charge center or toward a negative charge center. If one 

(text continued on page 64) 
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Figure 4.6 Simplified Drawing of the Expected Behavior of 
Triggered Lightning on the F106B for Ambient 
Electric Field Oriented Tail to Nose 
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assumes once again that the initial breakdown is at the nose and a possible 
second breakdown occurs at the aft end, the expected field behavior is 
shown in the figure. It looks very similar to the waveforms of Figure 4.3, 
simply reversed in polarity. The actual situation is somewhat different 
between the two cases however. In the case of Figure 4.3 a negative corona 
forms around the nose, since the front of the aircraft is negatively charged. 

In Figure 4.6 a positive corona forms there. It is well known that a negative 
corona requires a higher electric field to form than a positive corona, so 
the details of the two cases can be expected to differ somewhat. Intuitively 
it would seem reasonable that the smaller field required to initiate a po- 
sitive corona would result in a slower and smaller measured D-dot response 
as the corona expands outward from the aircraft. The negative corona would 
result in a faster breakdown when the required higher field is finally 
reached. This would appear to be characteristic of the measured data, in 
that initially positive D-dot records (near the F106B nose) tend to be larger 
and have more high frequency content than initially negative records. 

This type of simple analysis has been done for the geometry having 
electric field oriented along the fuselage, three possibilities of charge 
on the aircraft, and two breakdown locations. The three charge cases are 
+0 -0 and zero charge, where Q has a magnitude such that it can signi- 

ficantly alter breakdown characteristics and locations. That is, Q q cannot 
cause air breakdown by itself, but it can change the magnitude of the electric 
field which results in breakdown. The two possible breakdown locations are 
the nose and tail. Rather than presenting a picture of each case separately, 
the results of the analysis are presented in Table 4.1. For each case is in- 
dicated the direction of the electric field, the charge on the aircraft, the 
point of breakdown, the initial excursion of the D-dot and B-dot sensors, and 
the initial sign of the current at the breakdown point. The sign convention 
for the current was chosen so that a positive sign indicates current flowing 
off the aircraft. 

It is interesting to note that the majority of the measured records 
which include simultaneous D-dot and B-dot data show both being positive for 
the initial excursion. This occurs only three times in Table 4.1, and one 
of those three cases is unlikely. The two remaining cases both have a front 
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Table 4.1 


Parameterization Indicating Initial Sensor Signs 
for Differing Choices of Electric Field Orientation, 
Charge on Aircraft, and Air Breakdown Location 
(B-F - back to front, F-B = front to back, 

N = nose, T = tail ). 
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to back electric field, and breakdown at the nose of the aircraft. The net 
charges for these cases are zero and a negative charge. This is consistent 
with the expected negative charging of the aircraft from collisions with 
ice particles. 

Cases in which the ambient field is not oriented along the fuselage 
are far more difficult to analyze simply, because then aircraft resonances 
are much more important for the D-dot and B-dot responses. That is, it isn't 
clear what D-dot behavior will result from a wing to wing triggered lightning 
event. For that case the shape of the aircraft must be considered in detail. 

To examine possible differences between triggered and natural 
lightning responses of the F106B, consider Figure 4.7. That figure shows the 
approach of a negative leader toward the nose of the plane. The aircraft 
polarizes as in Figure 4.6 and the D-dot and B-dot responses as attachment 
occurs are also much the same as in that figure. Physically the only major 
difference between the triggered and natural lightning cases is the non- 
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Figure 4.7 Simplified Drawing of the Expected Behavior of the 

Fields on the F106B During the Approach and Attachment 
of a Natural Leader Channel 
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uniformity of the electric field around the tip of the leader. This would 
change the polarization of the aircraft before attachment somewhat, but 
the responses of the D-dot and B-dot sensors at attachment and detachment 
should not be greatly altered. Therefore, it seems unlikely that the two 

D-dot and B-dot records alone can distinguish between triggered and natural 
1 ightning. 

4.5 Effects of Triggered Lightning on Other Aircraft 

Because the NASA thunderstorm research program is eventually 
aimed at determining the effect of lightning on other aircraft than the 
F106B, it is of interest to examine whether the study of lightning trigger- 
ed by the F106B accomplishes that objective. Is triggered lightning, 
specifically that triggered by the F106B, of interest to other aircraft ? 

The answer to this is that it is of. interest, in that most lightning strikes 
to aircraft may be triggered. Although natural lightning strikes surely do 
occur, they are probably far less frequent, though possibly more severe, 
than triggered lightning. They will occur more frequently at low altitude 
than triggered strikes, because an aircraft at low altitude is unfavorably 
oriented to produce large local electric fields in the presence of a typical 
thunderstorm. Therefore at low altitudes a given strike has a greater prob- 
ability of being due to natural lightning, but one must also expect the 
probability of any strike occurring to be lower. 

Then given that the majority of strikes to aircraft could be from 
triggered lightning, certainly one must investigate its effects. The dif- 
ficulty is that triggering by a given aircraft is strongly dependent on 
the shape of that aircraft. For example, the F106B with its aerodynamic 
shape and sharply pointed nose is far more likely to trigger a lightning 
strike than is an aircraft of similar size designed for slower speeds. Also 
aircraft such as commercial jets and military transports, though not having 
sharp points, can trigger lightning simply by virtue of their large size. 
Hence each aircraft must be considered individually with regard to its shape 
to determine its potential for triggered lightning. 
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In addition, typical flight altitude must be considered, since the 
breakdown field is a function of height. At high altitudes two effects com- 
bine to increase the probability of triggered lightning. Because the air 
density is less the breakdown field is less, so the enhancement of the ambient 
field required to achieve breakdown is smaller. Secondly, fields at high 
altitude tend to be oriented more horizontally, and as such can be enhanced 
significantly by an aircraft in level flight, something which does not happen 
for the predominantly vertical fields at lower altitude. 

The advantage of triggered lightning over natural lightning is 
that it can be predicted to some extent. Enhancement factors for the field 
around a given aircraft can be easily calculated, and using these one can 
predict with a fair amount of certainty whether triggered lightning will 
occur in a given ambient field. This opens the possibility that it can be 
mitigated either by avoidance of the conditions that cause it, or by 
protecting the aircraft from its deleterious effects. 

There are other conditions not covered in the simple considerations 
of this chapter which can affect a triggered lightning event. One of these is 
the net charge on the aircraft. The net charge will always enhance the ef- 
fect of an ambient field at some locations on the aircraft and suppress it 
at other locations. Hence, depending on the sign of the charge and the 
direction of the ambient field, a net charge may either suppress or enhance 
a triggered lightning event. It can also move the location of the initial 
breakdown from one place to another. 

Another factor that has not been considered is the effect of 
particles in the air around the aircraft. These could be rain, cloud drop- 
lets, hail, graupel , or ice crystals. For most of the F106B data, ice 
crystals are most likely, because of the flight altitude and low air temper- 
atures. These particles, depending on conductivity, can locally enhance 
electric fields around themselves, and this could further support a triggered 
lightning event. This may be particularly important for the case of ice 
crystals, where irregular shapes are likely to be found. The effect of 
particles will be investigated in future work. 
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A third factor which has not been considered is the effect of air 
pressure variations around the F106B while in flight. Because air breakdown 
is strongly dependent on the local pressure, these pressure variations may 
be expected to significantly affect the location of initial air breakdown. 
These locations should show a strong preference toward low pressure regions. 
Also, because the pressure around the aircraft is a function of air speed, 
the presence or absence of triggered strikes is expected to be a function of 
flight speed. 
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CHAPTER 5 


TRIGGERED LIGHTNING MODELING 

In order to model triggered lightning it is necessary to place the 
F106B in a static electric field with an arbitrary orientation, and to allow 
for an arbitrary net charge on the aircraft. Physically what happens is that 
the F106B flies into a slowly increasing field and picks up a net charge over 
a period of seconds or minutes. This slow buildup means that at any given 
time the field and charge distributed on the aircraft are static. Because the 
triggered lightning model [4] operates with one nanosecond time steps, the 
slow buildup cannot be modeled. Therefore it is necessary to use the static 
solution as an initial condition. The way in which this initial condition is 
found is described below. 

The modeling procedure is as follows: Four linear codes are run 

initially. In three of these codes the ambient electric field is oriented 
along one of the coordinate axes, either along the fuselage, wing to wing, 
or vertically. Huygens' surfaces [6] are used to produce the field, which 
is brought from zero to one volt per meter in a time of four hundred nano- 
seconds. The long risetime of the field is needed to minimize the excitation 
of aircraft resonances. The code is then run to one microsecond at a constant 
one volt per meter ambient field level in order to assure that aircraft reso- 
nances have died away and a static solution is obtained. At the one microsecond 
point all the fields in the problem space are stored so that the electrostatic 
solution at that time may be used as initial conditions for the nonlinear 
corona model . 

The fourth code to be run introduces a known current which flows 
between the F106B and the boundary of the problem space. This produces a net 
charge which is then allowed to flow around the aircraft until equilibrium 
has been reached. At this time all the fields in the problem space are again 
stored to provide initial condi tons for the nonlinear model. This is then 
the linear static solution for a known net charge on the aircraft. 
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Because these four codes are linear the superposition principle 
applies, and any ambient field and net charge combination can be constructed 
from the results. Therefore, by combining the linear results properly, 
initial conditions for arbitrary orientation and amplitude of electric 
field and arbitrary charge on the F106B can be obtained for the nonlinear 
code. For example, consider the case of an ambient field of magnitude E 
oriented at a 45° angle between the left wing and the tail section of the 
F106B, with no vertical component. For simplicity assume there is no net 
charge on the aircraft. The 45° angle means that there are equal components 

of the field along the fuselage and in the wing-wing direction of magnitude 

’ 707 E o* The static solution for this case is constructed from two 
sol utions : 

1) solution with electric field along the fuselage; 

2) solution with electric field oriented from right to left wing. 

If one designates the first solution by S ] and the second by S , then the 
solution to the particular case being considered may be written, 

S = .707E q S 1 + .707E q S 2 . The solutions are first scaled from the 1 V/m 
normalized field to the desired magnitude, and then added together to get 
the solution to the more general case. In a similar manner, the static 

solution for any ambient field and net charge can be constructed from the 
four basic solutions. 

Finally, using these linearly derived initial conditions, the 
nonlinear code is run for a microsecond, which is normally long enough to 
observe air breakdown at one or more points on the aircraft. 


There is a point that should be made here about the correspondence 
between actual triggered lightning and the model. Before triggering, the 
assumption is made that the aircraft is flying through a slowly increasing 
ambient electric field, so that the field distribution can be regarded as 
static at all times. A slowly increasing field is one that changes very 
little in a typical air breakdown time scale (^few microseconds). This being 
the case, it is clear that the aircraft will form corona and probably trigger 
lightning at the lowest ambient field for which this is possible. That is. 
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given that the field is slowly changing, air breakdown will occur at some 
point on the aircraft as soon as the field at that point has reached minimum 
intensity for breakdown. No field on the plane can be above that minimum 
intensity, at least not before breakdown has occurred. After breakdown at 
some point has occurred, the field distribution may vary quite rapidly, al- 
lowing some fields to rise above minimum breakdown intensity. At any rate, 
the fact that it is always a minimum intensity field which initially causes 
air breakdown means that corona formation will occur quite slowly, possibly 
taking several microseconds to develop. Slow corona formation implies that 
triggered lightning D-dot and B-dot amplitudes measured on the F106B should 
be quite small when compared to those caused by the strike of a fully formed 
leader at some attach point. 

The relevance of this discussion to the nonlinear modeling is as 
follows. The nonlinear model uses as an initial condition what is essentially 
the linear solution to the problem of the aircraft in a static ambient field 
with an arbitrary charge on the aircraft surface. Because the solution is 
linear, fields can exist in the problem space which are above air breakdown 
intensity. This is not possible physically, because any field above breakdown 
level would of necessity have caused air breakdown at some earlier time while 
the ambient field was increasing. Hence the proper initial condition for the 
nonlinear model is a static field distribution in which the largest field is 
of air breakdown intensity. However, this is difficult to accomplish exactly, 
and it is often the case that at some point in the problem space the local 
field is larger than the minimum breakdown field. At that point corona will 
grow much faster than it would in reality, giving rise to predicted D-dot 
and B-dot responses that are much higher than are measured. The trick the 
analyst must perform is to reduce the ambient field until the largest field 
on the aircraft is at or near minimum breakdown intensity. For that case only 
can the predicted and measured responses be expected to agree. In the predicted 
waveforms to be presented later in this chapter, some amplitudes are far above 
any responses which were measured. That is because the ambient field chosen 
for those cases was too large, resulting in an unrealistic corona growth. 
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Linear enhancement factors for selected points on the block model 
of the F 1 06 B are shown in Figures 5.1 - 5.4. In Figures 5.1 to 5.3 the numbers 
are the electric field intensity in the direction of the arrow for an 

ambient field intensity of one volt per meter in the direction indicated on 
the figure. In Figure 5.4 the numbers are the electric field intensity in 
the direction of the arrow for a one microcoulomb charge on the aircraft. 

These numbers allow a rough estimation of how large the ambient field and 
charge need to be in order for an air breakdown field to exist at some point 
on the aircraft. Note that the largest field enhancement for the ambient 
field cases (Figures 5.2 - 5.3) is at the nose of the aircraft, as expected. 

In Figures 5.5 to 5.16 are presented the results of the triggered 
lightning runs which have been completed. They consist of D-dot, B-dot, and 
current waveforms as a function of time. The current waveform is the 
current flowing onto or off the aircraft at a location near the expected air 
breakdown point. For example, if the expected breakdown point is at the nose 
of the F106B, the current is that flowing on the fuselage just behind the nose. 
The environmental conditions for each run (e.g. orientation of field, charge 
on F106B, relative air density, humidity) are indicated in the legend of 
each figure. 

Several points should be made about the triggered lightning results. 
First, note that the largest ambient field runs produced predicted D-dot 
and B-dot responses which are much larger than any measured responses. The 
reason for this was discussed earlier in this chapter. Second, notice that 
when the ambient field level is reduced, the responses are slower and smaller 
in amplitude. This too is expected since the corona forms more slowly for the 
lower field levels. Although it is not apparent from the figures, the re- 
sponses also occurred later in time for the lower ambient field levels. Also, 
note that the predicted response amplitudes approach the measured amplitudes 
as the ambient field is reduced. Finally, one should be aware that the pre- 
dicted ambient field levels for triggering are certainly somewhat high. The 
block model of the F106B produces local fields near the sharp edges and points 
of the aircraft which are lower than would be found around the real aircraft 
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for the same ambient field level. Hence the ambient field needed to trigger 
lightning from, for example, the nose of the aircraft may be as much as a 
factor of two or three lower than that seen in Figures 5.5 to 5.16. This prob- 
lem will be investigated in more detail when the subgrid (Chapter 7) is in- 
corporated into the nonlinear model. 

Figures 5.17 and 5.18 show two of the predicted nonlinear responses 
from the same run overlaid on a pair of simultaneous D-dot and B-dot records. 

It is worth mentioning that the nonlinear triggered lightning model is the 
only model which has been able to even approximately match both the D-dot 
and B-dot records in amplitude and waveshape. With linear models it has 
been possible to match one or the other of the measured records exactly, 
but never both in the same run. This gives credence both to the hypothesis 
that triggered lightning is responsible for the majority of the measured 
F106B responses, and to the nonlinear corona model as a vehicle for analyzing 
the phenomena. 

The reader may note that all the currents calculated by the non- 
linear model are small, usually less than a kiloamp at peak. This is not only 
a function of the triggering condition (e.g. electric fields near breakdown 
intensity), but is also a function of the model. Because the current model 
cannot follow the formation of a lightning channel, the currents are restricted 
to levels appropriate to times before channel formation. It is hoped that the 
more complete model described in the next chapter will correct this shortcoming. 

In addition to the triggered lightning calculations for the F106B, 
one run has been done to investigate the triggered lightning response of the 
Cl 30 aircraft. Although the Cl 30 is larger than the F106B, it is considerably 
less streamlined, and therefore has fewer sharp points. Because of this the 
ambient field needed to trigger on the C130 is larger by almost a factor of 
two over that on the F106B. Figures 5.19 and 5.20 show the D-dot, B-dot and 
current responses on the Cl 30 for the case of ambient electric field oriented 
along the fuselage. Figures 5.21 and 5.22 show the corresponding frequency 
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Figure 5.17 Calculated Non-Linear D-Dot Overlaid on Measured D- 
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Figure 5.18 


Calculated B-Dot Overlaid on Measured B-Dot from 
Flight 82-042 (Ambient Field 1.9 x 10$ V/m, Oriented 
Nose to Tail) 
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transforms. The responses were calculated on the aircraft in approximately 
equivalent locations to those on the F106B. The biggest difference from the 
F106B responses is in the obvious resonant behavior of the C130. This is 
because the Cl 30 ' s lower frequency resonances are strongly excited by the air 
breakdown event. The widely variant behavior of the C130 and F106B responses 
illustrates the necessity of treating each aircraft individually when dealing 
with triggered lightning. Both the ambient field needed to trigger and the 
response after air breakdown are a strong function of aircraft geometry. 
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CHAPTER 6 


ENHANCED NONLINEAR MODELING 

This chapter reports on an extension to the nonlinear finite 
difference code which is currently being developed. The extension involves 
including fluid momentum and energy conservation equations to the particle 
conservation equations of the previous model. In that previous model 
heating of the electrons and ions in an air breakdown event was ignored, 
and organized motion of these charged particles was allowed for by an 
experimentally determined mobility. This mobility replaced the momentum 
conservation equation of the current model. Table 6.1 summarizes the ma.ior 
differences between the models. In what follows, the previous model will be 
referred to as model I and the newly developed model as model II. 


Table 6.1 

Differences Between 

Model I and Model II 

Physical 

Quantity 

Implementation in 
Previous Model (I) 

Implementation in 
Current Model (II) 

Charged particle 
number density 

Particle conser- 
vation equations 

Particle conser- 
vation equations 

Charged particle 
motion 

Mobility 

Momentum conser- 
vation equations 

Particle heating 

Ignored 

Energy conser- 
vation equations 


The question naturally arises as to what is gained by using 
model II rather than model I, especially in view of the fact that model II 
is significantly larger than model I in terms of computer requirements. 
There are many more physical quantities to keep track of, and the numerical 
integration of the momentum and energy equations increases computer run 
time considerably. The answer is that the newer model is more complete. 
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providing self-consistent information about drift velocities and particle 
temperatures. These new predictions make it easier to compare the model 
to experimental results, where particle temperatures are often measured. 

The newer model is also more versatile, since it provides a framework for 
still more complete models, in which more energy and momentum sources and 
sinks are considered. If desirable, one can add more species of particles 
to the model by differentiating between the different types of heavy par- 
ticles. This would allow one to implement more precise ionization, re- 
combination, and excitation rates for each species, rather than an average 
over the gas as a whole, as is done presently. 

Model II, like model I, is a fluid dynamic model, using fluid 
conservation equations to calculate charged particle densities, motions, 
and temperatures. The model contains three fluid equations for the electrons: 
conservation of particles, conservation of momentum, and conservation of 
energy. The heavy particles are a bit more complicated. There are particle 
conservation equations for positive and negative ions separately. There is 
a momentum conservation equation for the positive ions and an energy con- 
servation equation for all the heavy particles together. The reason there are 
no separate energy equations for the heavy particles is that these particles 
are very strongly coupled by collisions because of their similar masses, 
resulting in a very rapid equipartition of energy. 

The particle conservation equations for electrons, positive ions, 
and negative ions are shown below. At present the neutral gas particles 
are assumed to represent an unchanging background, so no equation is 
needed for them. 

(n eV = Q + 9 n e ‘ a e n e ' en e n + + k e 

(n + v + ) = Q + gn e - Bn g n + - Y n + n_ (6.1) 

(n_v_) = a g n e - Y n + n_ 


9t 

9n + 

3t 

3n_ 

w 


+ v • 


+ v 


+ v 
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These equations determine the time evolution of the charged 
particle number densities (n g , n + , n_) at a given point in space. The left 
hand side of each equation is the total time rate of change of the density 
at a point and takes into account organized motion of the particles through 
the convective derivative. The right hand sides represent sources and sinks 
of particles. Q is the ambient ionization rate from such sources as cosmic 
rays. gn e is the ionization rate due to avalanching in the electric field. 
This provides electrons and positive ions and depends on the strength of 
the electric field and the air density. a e n g represents attachment of 
electrons to neutral particles to form negative ions. Hence it represents 
a sink for electrons and a source for negative ions. 6n g n + is the rate of 
electron positive ion recombination, and ^n + n_ is negative ion-positive 
ion recombination. The last term in the electron density equation, k g , 
represents diffusion of electrons from one location to another because of 
a temperature or density gradient. 

The details of most of these terms have been documented in 
previous reports [3,4]. k e has not been documented before, so a simple 
derivation of its form will be given here. Consider a surface separating 
two volumes of fluid with different densities and temperatures as shown 
in Figure 6.1. 





Figure 6.1 Two Equal Volumes of Fluid Separated by a 

Mathematical Surface for Use in Calculating 
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It will be assumed for simplicity that all particles in a given volume are 
moving at the speed appropriate to that temperature. In reality, of course 
there is a distribution of velocities, but this merely complicates the 
analysis, without contributing to a physical understanding of the process. 

The first step in the analysis is to calculate the net flux of 
particles across the surface. Only the velocities in the x direction need 

be considered, since the other components do not contribute to the flux. 
Hence, 


1 2 i 

2 rnv x-| " 2 kT 1 

1 2 i 

2 mv x 2 = 2 kT 2’ 


( 6 . 2 ) 


where m is the particle mass and k is Boltzman's constant 

v 

x 1 V m 



(6.3) 


To calculate the flux it must be assumed that half of the particles in 
each volume are moving in the +x direction and half in the -x direction. 
Then the number of particles crossing the surface per unit area per unit 
time can be written. 


F = 2 n l V X] ' 2 n 2 v x ? = ("l A " n 2 A )• (6 - 4 ) 

The change in density per unit time because of this flux can now be 
written. 


An _ r aA _ , i IT . / 

At- F AV = ± 2 W ("lA^A)’ (6 ‘ 5) 

where aA is the area of the surface, and it has been assumed that V = V 
= AV. The ambiguity in sign refers to the fact that one of the volumes 2 
gains density while the other loses density. In differential form Equation 
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(6.5) can be written. 



diffusion 



3 _ 

9X 


(n /f). 


( 6 . 6 ) 


This is the desired form for k g . Note that it depends both on a density 
and temperature gradient. The effects of each can be separated by expanding 

the derivative, 

(6.7) 


"| /k / x n \ 

k e = - 2\lm ' + 2 /f 3X ' 


Here the tirst term in parentheses is the effect of a density gradient 
with uniform temperature, and the second is the effect of a temperature 
gradient with uniform density. Note that k g depends on nf , and is 
therefore much more important for electrons than for the heavy particles. 
That is why Equations (6.1) only include a diffusion term for electrons. 

The momentum conservation equations for the three species are 


shown below* 


n 
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3V 
a 

dt 




n q 

a a 

m 




— vp 
m a 
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v 

a 


o 

C 


( 6 . 8 ) 


Here a refers to species type; electrons or positive ions. 

Model II ignores the motion of negative ions, making the assump- 
tion that virtually all of the transfer of negative charge is accomplished 
by electrons. This is especially true during the active early portions of 
a discharge. Only when the discharge is dying away does negative ion 
motion become important. 

Note that Equation (6.8) is a vector equation, and that the con- 
vective derivative on the left hand side couples the different components 
of v together. This term is an additional source of nonlinearity in the 
problem, beyond the introduction of electric field dependent conductivity. 

The left hand side of Equation (6.8) is the total rate of change 
of the fluid momentum per unit mass, including that due to organized par- 
ticle flow. The first term on the right is the force per unit mass exerted 
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on the fluid by the electric and magnetic fields. Because of the factor 
v q /c, where c is the speed of light in a vacuum, the magnetic force term 
is expected to be much smaller than the electric force term. The second 
term on the right represents the force per unit mass exerted on the fluid 
by the partial pressure gradient. The force is directed from high pressure 
regions to low pressure regions, thus tending to equalize the pressure. 

The third term on the right is the force per unit mass exerted on the 
fluid because of collisions. At the present time only collisions with 
neutral particles are considered because of the low level of ionization, 
and the collision frequency is assumed to be a constant. These approxi-’ 
niations can be improved in the future if results warrant. The force due 
to collisions is always oriented opposite to the direction of the velocity. 
Hence the drift velocity in steady state) is usually a balance between 
the driving force of the electric field and the damping of the collisions. 

A crude approximation to v a drift is found to be. 


drift 

a 


m u 
a C 


t . 


(6.9) 


Note that in model I a mobility was defined such that v dnft = u t Then 

a a 

Equation (6.9) implies that an approximation to u q is . This approxi 
mation can also help to fix u if u is known ° 

C a 

The energy conservation equations are the most complex of the 
fluid equations used in model II. At present there are only two, one for 
electrons and one for heavy particles, but it may be necessary to include 
more eventually. Energy reservoirs that may need to be considered are 
vibrational and rotational states of molecules, and excited states of 
molecules, atoms, and ions. These states can take energy from particle 
motion, retain it for some period of time, and then release it as kinetic 

or electromagnetic energy. The present model II only begins to account 
for these effects. 


The present energy balance equations are shown in Equations 
(6.10). As was stated earlier, the reason the heavy particles are all 
lumped together in one equation is that they are strongly coupled by 
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collisions, resulting in rapid equipartition of energy. In an electron- 
heavy particle collision, the electron retains most of its energy, hence 
electrons may have a different temperature than the heavy particles. 
However, when two heavy particles collide, a large fraction of the kinetic 
energy available can be transferred. Therefore, equipartition of energy 
occurs in at most a few collision times. This is much less than the time 
step used for model II, so the heavy particles can all be considered to 
have the same temperature, and one energy balance equation used for all. 


£ * (v e -v) c e - Ve t . % - J u c (., - S,°> * 9"e (7 m e V 2 - E ion ) 


* e q ' “ e E e ' 6n * E e + H 


e e excitation 


(6.10) 


rn 


^ -V) £ H = <4 K + "-) 1 + ST u c < 


Ox 

£ e - e e } 


1 2 

- 2 9 m e n e v + + ^excitation" 


Each of the terms in the above two equations will now be discussed indi 
vidually, with the approximations used in each indicated. 

The left hand side of each equation is the total rate of change 
of the energy density of the species. This includes changes due to organized 
flow of the particles, through the convective derivative term. The energy 
density is defined as the sum of the energy densities included in organized 
motion (particle flow) and internal energy (temperature). Note that the 
heavy particle equation uses the velocity of the positive ions, v + , in the 
convective derivative and throughout the rest of the equation. This is done 
because the positive ions appear in the early stages of air breakdown, 
before the electrons are able to attach to neutrals to form the negative 
ions. Hence it is necessary to use v + in the first part of a discharge, 
and the approximation is made that this is the flow velocity of the heavy 

fluid at all times. 

The first term on the right in the electron equation, q g n e t • v g 
is the energy density transferred to the electron fluid per-unit time by 
the action of the electric field on the electrons. Because v e is nearly 
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always oppositely directed to £, and q g is negative, it is generally the case 
that energy flow is from the field to the electrons. However, it is at least 
mathematically possible (with large pressure gradients, for instance) for 
Vg to be in the same direction as It, in which case energy is transferred from 
the electron fluid to the field. It may be noted here that the magnetic field 
makes no contribution to the energy density because of the nature of the force. 
The energy transferred from the magnetic field contains a term of the form 
(v e x I) • v e which is, of course, identically zero. 

m 

The second term on the right of the electron equation, — u (e -c °) 

1 m c e e ' 

represents energy transfer from the electrons to the heavy fluid H 
because of collisions. These are electron-neutral collisions in general since 
the degree of ionization is assumed to be quite low, at least initially. In 
this term u c is the collision frequency, and e e ° is the ambient electron 
energy density. Eg 0 is a lower limit for the electron energy density, below 
which is not allowed to drop. Note that the energy transfer due to col- 
lisions is proportional to the electron energy, the collision frequency, and 
the ratio of the masses. The latter proportionality comes about because of 
the dynamics of a single electron-heavy particle collision. Less energy is 
transferred as the ratio of the masses becomes smaller. An approximation is 
made in this term that all heavy particles can be assigned some average mass, 
m H . 

The third term on the right, gn g (1 m e v + 2 - e iQn ), represents the 

energy gained (or lost) by the electron fluid because of avalanching. In 

avalanching, electrons are liberated from neutral particles by collisions 

between the neutrals and energetic electrons. The energetic electron loses 

an amount of energy equal to the ionization potential of the neutral, e. . 

... ’ion’ 

and it is assumed that the new electron is born with a kinetic energy 

appropriate to the motion of the heavy fluid. Both of these energies are 

approximations, as in reality, a statistical distribution is to be expected 

in the energetics of ionizing collisions. 

Eq represents the energy gained by the electron fluid because of 
the ambient ionization rate. This may be an energy over and above that added 
simply by the liberation of an electron. For instance, cosmic ray ionization 
may result in an electron released with a large amount of kinetic energy, 
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which is eventually distributed throughout the electron and heavy particle 
fluids. 

The next two terms in the electron equation, a g e e and en + e g , 
represent the energy lost to the electrons because of attachment to neutral 
particles and because of electron-positive ion recombination. Both of these 
processes are assumed here to be a simple loss of the total energy of the 
electron engaged in the process. That is, for example, in the case of an 
attachment, the energy of the electron involved is assumed to be released 
in the form of a photon, which is then lost from the gas. That is not 
necessarily the case, because the energy may be transferred to heavy par- 
ticle kinetic energy or stored in an excited state of the heavy particle. 

The probability of these processes may be explored in future enhancements 
to model II. 

The term H g in the electron equation is not a source or sink term 
as such but represents a transfer of energy from location to location ana- 
logous to the k g term in the particle conservation equations. Its form is 
similar to the form derived there. The term is nonzero only in the presence 
of temperature or density gradients. 

k .. . represents a process which has not yet been included 

exc i xs l i on 

in model II, but is felt to be essential to its proper functioning, and will 
be added shortly. The process is excitation of vibrational modes of molecules 
in the gas by electron collisions. The vibrational modes eventually decay 
into kinetic energy of the molecules, hence providing a significant path 
for energy transfer from electron to heavy particle kinetic energy. Direct 
collisions cannot accomplish this because of the great disparity in mass 
of the two types of particles. In results to be presented later it is seen 
that calculated electron temperatures are quite a lot higher than is seen 
experimentally. This occurs because in the present model II the electron 
fluid has no efficient way to transfer its internal energy to the heavy 
fluid. The addition of k exc -j tation t0 the model wi11 remedy this and hold 
down the electron temperature. The actual form to be used for k exc i tat i on 
is yet to be developed. Initially, at least, it is expected 
that the transfer ot energy through the vibrational modes will be in- 
stantaneous, with no energy stored in the modes. This may be changed to 
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allow some time delay from excitation of the modes to deexcitation if this 
is found to be necessary. 

On the right hand side of the heavy particle energy equation, 
it can be seen that many of the same terms as in the electron equation 
are present. The only new term is the first on the right hand side, 

q + (n + + n_) t ■ v + . This represents the energy provided to the heavy par- 
ticles by interaction with the electric field. Notice that this includes the 
contributions from both positive and negative ions. Strictly speaking, the 
contribution from the negative ions should be written q n. £ • v , but if 
l v _l % |v + |, and q + = -q_, the form shown in Equations (6.10) is very nearly 
correct. 

Equations (6.10) are a simplification of the actual energy balance 
in an electrical discharge. As stated earlier, there are many other processes 
and energy reservoirs that could be included in model II. Some of these in- 
clude energy stored in vibrational molecular states, excited atomic and 
molecular states, and a photon fluid arising from the statistical decay of 
excited states and recombinations. These have not been modeled at the present 
time in order to keep the model as simple as possible. They can be added in 
the future if it becomes necessary. 

Figure 6.2 shows the physical situation to which model II has been 
applied. This is similar to the experimental setup of Collins and Meek [7 ] 
which has been investigated in previous reports [3,4]. However, the apparatus 
of Collins and Meek was significantly larger in the radial dimension. A 
smaller problem space has been used here because of run time considerations. 

The problem space shown in Figure 6.2 was gridded with a spatial 
resolution of one centimeter in both the radial and axial directions. The 
time resolution was 20 picoseconds. The source was a voltage of 7 x 10 
volts applied across the plates having the form of a step function with a 
risetime of 100 nanoseconds. Air parameters chosen are shown in Table 6.2. 
These are intended to be representative rather than exact values. For 
example, the mass of the heavy particles in air is not well defined, so the 
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Figure 6.2 Schematic of Problem Space Modeled by Enhanced 
Nonlinear Air Chemistry Code 


Table 6.2 Parameter Values Chosen 
for Model II Results 

Relative air density = 1. 

Percentage water content = 6. 

1 3 -1 

Collision frequency = 1 x 10 sec . 

“26 

Heavy particle mass = 2.34 x 10 kg. 

“18 

Ionization energy = 2.4 x 10* joules. 

Ambient temperature = 300°k. 

24 -3 

Ambient neutral particle density =1x10 m . 


108 



mass chosen is intended to be in some sense an average value. 


The results from model II are shown in Figures 6.3 to 6.9. In 
Figure 6.3 is presented the electric field behavior at the tip of the rod. 
The reader will note that the field at breakdown drops essentially to zero 
and then rises again at a later time. The second rise occurs because of the 
constant voltage source across the plates. After the initial breakdown, the 
electrons which provide the local conductivity begin to attach to neutrals 
to form negative ions and also to recombine with positive ions. This drasti- 
cally lowers the conductivity and allows the voltage source to put energy 
back into the electric field. 

Figures 6.4 and 6.5 show the radial and axial electron velocities, 
respectively, near the tip of the rod. Note that these drop to zero as the 
electric field disappears. Figure 6.6 shows the electron temperature. This 
figure is somewhat deceptive, as the temperature is large at times when the 
electron density is small. The temperature is not really a significant 
quantity when there are only a few electrons present. Nevertheless, the 
predicted electron temperature is too high, and it is expected that the 

of the kg XC -j tg.j-.jQp term in the electron and heavy particle energy 
equations will hold the value within acceptable limits. 

Figure 6.7 shows the electron density at the tip of the rod. 

Notice the growth due to avalanching and the almost equally rapid decay 
because of attachment. The positive ion density is shown in Figure 6.8. Here 
again the growth is because of avalanching. However, the decay is from 
recombinations, a much slower process than attachment, so the positive ion 
density stays high. Figure 6.9 shows the heavy particle temperature, which 
is virtually unchanged from the ambient value. This is another indication 

that Excitation is needed in the energy equations to transfer energy from 
the electrons to the heavy particles. 

Future plans for model II include a study to determine which terms 
are most important in the conservation equations. This will possibly allow 
a reduction in size of the final version. The model must also be checked in 
as many experimental situations as possible to increase confidence in its 
validity. Finally, it would be desirable if model II could in some way be 
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FIGURE 6 5 THE ELECTRON FLUID AXIAL VELOCITY NEAR THE TIP OF 
THE ROD 












incorporated into the three dimensional code containing the F106B. There 
are two possibilities for accomplishing this. The first depends on the 
final size of model II. If it is small enough, the model can be made three 
dimensional and incorporated directly into the nonlinear finite code. In 
this case it may be that only selected regions (nose, wing tips, etc.) 
would be considered for nonlinearities. The second possibility concerns 
interfacing the present two dimensional model II with the three dimensional 
code. This could only be done around regions of the F106B which exhibit 
two dimensional symmetry, such as the nose. Details of how the interface 
would work are not yet available. 
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CHAPTER 7 


SUBGRID THEORY AND DEVELOPMENT 
7.1 Introduction 

In standard applications of explicit finite difference techniques 
to solving Maxwell's equations the fundamental limit of spatial resolution 
is the largest spatial grid interval. In fact, because of problems with 
error propagation at short wavelengths, the rule of thumb limit of spatial 
resolution is normally taken to be that approximately five of the largest 
spatial grid intervals are required to resolve the wavelength of the 
highest frequency of interest. While this resolution is usually sufficient 
for problems involving the large scale response of a scattering object to 
electromagnetic fields, it is known that such coarse resolution results 
underestimate the magnitude of the charge densities and currents on a 
scatterer in the vicinity of sharp edges or points such as the nose of an 
aircraft. One way to understand this is to consider the finite difference 
solution as a volume average over the coarse grid unit cell. The larger 
this volume is, the lower the volume average is relative to the continuum 
results near a discontinuity. 

This is of particular importance when applied to the nonlinear 
corona model of the F106B. Corona and air breakdown are strongly dependent 
on the absolute magnitude of the electric field in a local region. Because 
the pitot boom of the F106B is very sharply pointed, large electric fields 
will exist near the tip. These large fields, however, will not extend very 
far from the tip. In fact the field will vary approximately as 1/r , where 
r is the distance from the tip of the boom. This is a very rapid variation 
which is impossible to model accurately in a nonlinear code with cell 
sizes of .5m x lm x .5m. This is not to say that air breadown will not 
occur in this model. It will, but it will occur at a larger ambient field 
level for the model than for the real aircraft. This may alter the D-dot 
and B-dot waveforms as predicted by the code. Therefore, it is necessary 
to attempt to achieve a better spatial resolution of the sharp points and 
edges of the aircraft. This needs to be accomplished without increasing 
the resolution at all points, because that leads to codes which become 
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rapidly unwieldy in size. One solution to this problem is to grid a portion 
of the aircraft (such as the nose) with a fine grid and then integrate that 
fine grid into the rest of the code. In many cases the nonlinear portion 
of the code can fit entirely into the fine grid, with the larger grid being 
linear. The details of this "subgrid" model are explained in the rest of 
this chapter. 

7.2 General Description 

In the model discussed in this report it is necessary to consider 
a rectilinear volume of space within which one wishes to obtain a finite 
difference approximation (f.d.a) to the solution of Maxwell’s equations. 
Well posed solutions require the specification of the vector fields over 
the entire volume initially and specification of the boundary values of 
the fields at all times. The f.d.a. used is an explicit, second order, 
centered difference procedure. In this procedure the Cartesian components 
of the vector and JH fields are assigned staggered grid locations. The 
bounding planes of the problem volume, are assumed to coincide with sur- 
faces containing the H field components tangential to the bounding surface. 
Contained within this volume is a small rectilinear volume in which the 
finite difference grid intervals are some even integer ratio of the larger 
problem volume grid interval. The bounding surfaces of the smaller volume 
also coincide with surfaces containing both the coarse grid resolution and 
the fine grid resolution tangential H field components. Conceptually, the 
boundary between the two volumes is the same as the inner boundary of the 
larger problem volume. The fields on this boundary must be specified at 
all times in order to obtain a well posed solution to either the coarse or 
fine resolution f.d.a. 's to Maxwell's equations. The details of these 
specifications are discussed in the section on boundary conditions. The 
essential point of the specification is that it must also be consistent 
with both the coarse and fine grid solutions. 

7.3 Details of the Fine Resolution Field Lattice and its Relationship 

to the Coarse Resolution Field Lattice 

As mentioned in the previous subsection, the boundary between the 
coarse and fine resolution grids is equivalent to any other boundary of the 
problem volume. In effect, it is the interface between two separate f.d.a. 
solutions to Maxwell's equations. 
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In the volume containing the fine resolution grid, the f.d.a. im- 
posed is again a centered difference staggered grid lattice. The grid inter- 
val is the fine resolution grid size. In Figure 7.1, a diagram of this 
lattice for a single field component (E z ) is shown in relation to a coarse 
grid field component. The expansion factor is 2, with the origin of the fine 
resolution grid corresponding to the origin of a given unit cell in the 
coarse grid. For this expansion factor, the x-y planes (planes normal to 
the z axis) between two full coarse grid x-y planes contain fine grid fields 
which are shared by coarse grid fields above and below this plane. Only those 
fine grid fields on the same x-y plane as the coarse grid are associated with 
a single field in the coarse grid. 

In order to avoid the uncertainties involved in extrapolating the 
boundary fields in time, both the coarse grid and fine grid fields are eval- 
uated at the same times. The Courant stability criterion then requires a 
time step which is stable for the fine resolution grid. 

It has been found that to obtain equivalent coarse grid fields from 
their fine grid counterparts, it is necessary to do a surface average. 

Consider a particular coarse grid electric field which is within the subgrid 
and therefore also represented by fine grid electric fields. For concreteness, 
assume that it is oriented in the X direction. The appropriate method for 
deriving the coarse grid field from the fine grid fields associated with it 

is the following: 

1) Determine the y-z plane which is at the same X location 
as the coarse grid field. 

2) Arithmetically average those fine grid fields which are 
located on that same plane and within a half cell spacing 
of the coarse grid field in the y and z directions. 

This average is the appropriate value to use in the coarse grid finite 
difference scheme. In equation form, the average can be expressed, 

F coarse _ / ] \ y f fine assoc -j ate d with the coarse grid field (7.1) 

x XX, ' *-• x 
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where and X z are the expansion factors in the Y and Z directions 

and p^ coarse an( j f^' ine are general coarse and fine grid x directed fields, 
respectively. The effect of this averaging process is to reduce the 
amplitude of Fourier components with wavelengths smaller than the coarse 
grid interval in coarse grid fields. 
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Figure 7.1 Distribution of Coarse (Z) and Fine (z) Grid Ez Fields for 
an Expansion Factor of 2 in all Three Coordinate Directions. 
The Solid Cube is a Coarse Grid Unit Cell. The Four Dashed 
Horizontal Planes Represent the Vertical Plane Boundaries 
of the Fine Grid Unit Cells. The Crossed Dashed Lines on 
Each Plane Correspond to the Projections of the Other Fine 
Grid Unit Cell Boundaries 
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7.4 Details of the Specification of H Fields Tangential to the 

Fine Grid - Coarse Grid Interface 

The general philosophy for specification of FI fields tangential 
to the coarse grid-fine grid interface is as follows: 

(i) Obtain the coarse grid H fields (tangential to the inter- 
face) at time t = n a t by using the coarse grid advance 
equations . 

(ii) Interpolate the coarse grid FI fields to obtain the 
fine grid tangential FI fields on the interface. 

(iii) Advance the fine grid H and E fields in the interior. 

(iv) Volume average the fine grid E fields tangential to the 
interface to obtain corresponding coarse grid E fields. 

This procedure completes one time step. The coarse grid determines 
the boundary for the fine grid by providing a distribution of tangential H 
fields which are used to obtain the fine grid FI fields by interpolation. The 
fine grid then couples back into the coarse grid by providing volume averaged 
E fields which are used in the coarse grid H field advance equations on the 
next time step. 

It is necessary to distinguish between interfaces which are penetrated 
by a scatterer and those which are not. When a scattering object crosses the 
fine grid boundary special care must be exercised in determining the H fields 
tangential to the interface. 

For interfaces which are not penetrated by a scattering object it 
is usually sufficient to perform a simple four point planar interpolation to 
obtain the components tangential to the surface. In the current implementation 
of the interface interpolation procedure the four point interpolation used 


is: 

f fine (x,y) = D + A < x -* 0 ) + B (y-y 0 > * c (*-* 0 >(y-y 0 ) 

(7.2) 

where 

*o * x< V^cfcoarse) 8 y o £y<y o + 4y c 


and 

D ■ F c<V y o> 
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It is clear by implication that these boundaries assume that the coarse grid 
field variation is relatively smooth (i.e., bilinear), and that the fine 
grid will follow this behavior. This is not the case when a scatterer penetrates 
the boundary. In order to match solutions across the coarse-fine grid inter- 
face, it is necessary to force a linear variation in fine grid magnetic fields 
inward into the fine grid subvolume. This is done by a linear interpolation 
of the outermost subgrid magnetic fields between the boundary and the corre- 
sponding subgrid field which is an expansion factor number of cells inside 
the subgrid volume. 

When a scatterer penetrates the interface the four point inter- 
polation procedure does not accurately predict the tangential H field's mag- 
nitude or distribution near the surface of the conductor in the fine resolution 
grid. In order to obtain a better estimate of the tangential H fields, a 
method that uses Maxwell's equations is required. For cases where the field 
distribution is approximately two-dimensional, that is 



one makes use of Euler potentials to obtain a solution for the distribution 
of tangential H fields in the fine grid at the interface. 

Define two scalar fields a(x,y) and P(z) = z- 

then H = vx( a ve) = -z x 7a = |y x - ~ 9 (7.4) 
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also 


2 • ( v x H) 


a 2 a 


ax 


ay 


aE 

2 

at 


( 7 . 5 ) 


or 


-2 J 

3 a , a a 

2 2 
ax ay^ 


- e 


5. 

at 


( 7 . 6 ) 


By specifying Neumann boundary conditions sufficiently distant from 
the surface of the conductor and interpolating to obtain the source terms for 
the fine grid resolution the c's can be solved for at the normal E-field 
positions in the fine grid. Note that the value of „ on a conductor is set 

to zero. At present, the Poisson equation is solved by the method of success- 
ive over-relaxation. [8 ]. 

Given the boundary values of the tangential H fields in the fine 
grid the interior points in the fine grid are advanced. The coarse grid E 
lelds contained entirely within the fine grid are then computed by taking 
the volume average of their associated fine grid fields. These fields are 

then used to advance the coarse grid H fields tangential to the coarse grid - 
fine grid interface. 

Gridding of Scatterers in the Fine Resolution Mesh 

When gridding a scatterer in the fine resolution grid some care 
must be exercised, particularly at the interface between the coarse grid and 
ine grid. The essential point is that all fine grid fields which contribute 
to the volume average of a coarse grid field must be zeroed. Hence all fine 
grid fields contained within a unit cell of a zeroed coarse grid E field 
as well as those shared by a zeroed coarse grid field and one which is not’ 
zeroed must be zeroed. It is often useful to diagram the fine grid and coarse 

grid fields in order to determine which fields are zeroed, since complicated 
scatterers may be hard to automate. 


7.5 
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7.6 


Resul ts 


An example of the utility and efficacy of this technique is a 
comparison of the normal electric field above one end of a 9m x 3m x 3m bar 
in a pure coarse grid and the same fields for a hybrid coarse-fine grid with 
the bar tapered to a lm x lm cross section. The pure coarse grid has a 
spatial grid interval of lm in all three Cartesian coordinates. The coarse 
resolution part of the hybrid has the same grid interval while the fine 
resolution part has a spatial grid interval of .5m in all three Cartesian 
coordinates. In all cases the bar is illuminated by a normally incident 
electromagnetic plane wave which rises to a constant 50 kV/m in 40 ns with 
a sin 2 leading edge. The incident electric field is polarized parallel to 

the bar's length. 


Figure 7.2a shows the waveform associated with the normal electric 
field in the pure coarse grid solution, whiie Figure 7.2b shows the normal 
electric field for the hybrid coarse-fine grid solution. Figure 7.3 
the waveform for the normal electric field when the bar cross section ,s 
tapered to lm. Note that the hybrid enhances the normal electric field 
cause of the improved resolution. The tapered bar could not be gndde 
in the pure coarse grid, but this is possible in the hybrid. 


The subgrid model has not yet been incorporated into the F106 
finite difference model or corona model. This will occur during the next 

year's effort. 
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Figure 7.2a Normal Electric Field Above Figure 7.2b Normal Electric Field Above 
the End of the Bar in the the End of the Bar in the 

Pure Coarse Grid Solution Hybrid Coare-Fine Grid 

Solution 



Figure 7.3 Normal Electric Field Above the Tapered 
End of the Bar in the Hybrid Coarse-Fine 
Grid Solution 
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CHAPTER 8 


F106B TRANSFER FUNCTIONS FOR THREE CURRENT 
INJECTION CONFIGURATIONS 

As an aid to a possible current injection test of the F106B, 
transfer functions have been calculated for three different configurations 
of the aircraft. The transfer functions were calculated between a current 
injected at the nose of the F106B and the D-dot and B-dot responses at the 
positions of the sensors. The three configurations used were: 


i) aircraft in free space, 

ii) aircraft resting on perfect ground plane, 

iii) aircraft elevated ten meters above a perfect ground plane. 

The purpose of the proposed test was to determine the transfer 
function between an injected nose current and the 0-dot forward and B-dot 
longitudinal responses of the F106B in flight. This of course, was not pos- 
sible experimentally, so the calculations were performed in order to determine 
what the effect the presence of the ground plane had on the transfer function, 
and whether elevating the aircraft ten meters above the ground significantly 
removed the effect of the ground plane. Also, in principle, the calculations 
allow the free space transfer function to be determined from the measured 
transfer function in the presence of the ground plane, because the differences 


are known. 

The D-dot and B-dot responses were calculated using a three- 
dimensional finite difference code including a thin wire channel for current 
injection [6]. The wire injected current at the nose of the aircraft and 
removed it from the base of the tail. The three configurations of aircraft, 
thin wire, and ground plane are shown in Figure 8.1. In the free space 
configuration, the thin wire runs directly into the problem space boundary 
and is driven by a voltage source at the edge of the problem space. The 
intent of this is to make the thin wire appear to be infinitely long on 
both the front and back of the aircraft. For the other two configurations 
the thin wire was connected to the ground plane on both the front and back 
in order to use the ground as a return path for the current. The voltage 


128 



Boundary of 
Problem Space 



(a) Free Space 



(b) Resting on Ground Plane 
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(c) Ten Meters Above Ground 
Plane 


Figure 8.1 The Three Current Injection Configurations 

Used in the Calculation of Transfer Functions 
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source for these configurations was placed at the junction of the wire and 
the ground at the front of the aircraft. In all cases the thin wire had^a 
characteristic impedance of 96ft and a velocity of propagation of 3 x 10 
m/sec. 

The transfer functions were calculated between the D-dot and 
B-dot responses and the current injected at the nose of the aircraft accord- 
ing to T = R/I, where T is the transfer function, ft is the response in 
question, Tf is the injected current, and tildes denote Fourier transforms. 

In the following the subscripts FS, GP and EL refer to the free space, 
ground plane, and elevated above ground plane cases, respectively. The 
transfer functions then may be designated T ps (D-dot), T ep (D-dot), T EL (D-dot), 
T FS (B-dot), Tgp(B-dot) , and T EL (B-dot). These transfer functions are, of 
course, complex functions of frequency. The frequency dependence is suppressed 
in the notation. 

Although the transfer functions themselves were calculated, it is 
more instructive to present ratios of the transfer functions. The ratio of 
the GP and EL cases to the FS case graphically illustrates the distortion 
in the transfer function caused by the ground plane. Figures 8.2 - 8.5 show 
the four possible ratios of the magnitudes of the transfer functions. The 
vertical scale is in dB and the horizontal scale is the base ten log of the 
frequency from 100 kHz to 50 MHz. Note that this means that the value of 
zero on the graphs indicates that the transfer function magnitudes are 
equal at that frequency. The figures also show that, as expected, the 
ratios are larger for the GP case than the EL case by amounts ranging up 
to 20 dB. Therefore, the effect of the ground plane is much reduced by 
elevating the aircraft, but it should also be noted that ratios of the 
order of 12 dB still exist for the EL case. The largest peaks in the ratio 
graphs are at approximately 2 MHz for Figures 8.2 and 8.3 and at 1 MHz 
for Figures 8.4 and 8.5. These appear to be characteristic frequencies of 
aircraft, wire, ground plane systems for the two configurations. As a 
caution in using the graphs, it must be remembered that the presence of 
a dielectric beneath the aircraft, such as wood support structure, could 
alter the transfer functions. The calculations were performed assuming the 
support had the same electrical properties as air. 
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CHAPTER 9 


SUMMARY AND CONCLUSIONS 

The work reported here has concentrated on analysis of the 1982 
inflight lightning data and the development of analytical tools to better 
understand that data. Data classification has been made in the time domain, 
involving the sorting of waveforms into groups having similar characteristics. 
Although somewhat subjective, this process may give some insight into similar 
lightning events that could be occurring on the F106B. Data analysis has also 
been done in the frequency domain by a numerical search for recurring reso- 
nances in the measured responses. 

The possibility that triggered lightning is responsible for much, 
if not all, of the direct strike data has prompted an investigation into this 
area. The expected behavior of triggered lightning and its effect on the 
F106B have been studied both intuitively and numerically. Calculated responses 
have been compared with actual measured responses. 

A pair of new analytical tools have been developed to help in the 
data interpretation process. The first is an enhanced nonlinear model, which, 
when finalized, will predict in much greater detail the evolution of a static 
discharge. Also developed was a subgrid model, which allows the analyst to 
spatially resolve a portion of a finite difference problem space in much 
greater detail than was previously possible. It is expected that both of 
these tools will be further developed and applied in following efforts. 
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APPENDIX A 


D-Dot and B-Dot Records 


The legend beneath each plot gives the flight and run number, 
strike attach and detach points, strike altitude, and ambient temperature. 
The attach and detach points were determined by NASA and LTI personnel 
from pilot's comments and the appearance of pit marks on the aircraft. 

The code for attach-detach locations is the following: 

U - Unknown 
N - Nose 
T - Tail 
LW - Left Wing 
RW - Right Wing 
M - Multiple 
N - Nearby 
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